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I.  INTRODUCTION 


Several  topics  under  the  general  title  of  the  progress  report 
were  investigated  during  the  time  period  from-  October  1,  1979  to 
September  31,  1980.  Essentially  five  topics  have  been  studied  during 
this  time  period;  three  are  essentially  experimental  investigations,  and 
the  other  two  are  theoretical.  The  five  topics  include: 

1)  Propagation  of  microwaves  along  a  plasma  column  and 
harmonic  generation  of  electrostatic  ion  cyclotron 
waves , 

2)  Imploding  tube  experiment, 

3)  Low  frequency  CW  (Hz)  rf  (IMHz-lOMHz)  driven  dc 
plasma  current, 

4)  Wave-particle  interaction  at  cyclotron  resonance,  and 

5)  turbulent  interaction  between  waves  and  charged 
particles. 

Topics  l)-3)  comprise  the  experimental  investigations  while  the 
last  two  topics  are  theoretical  in  nature.  In  the  next  section  we 
discuss  the  progress  to  date  on  each  of  these  topics. 
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II.  PROGRESS  REPORT  AND  WORK  TO  BE  CONTINUED 


1)  Imploding  Tube  Experiment:  Project  is  finished.  A  tech¬ 
nical  report  is  issued:  (Charles  D.  Hechtman  and  B.R.  Cheo,  "Inves¬ 
tigation  of  Electrically  Exploding  Thin  Shelled  Metallic  Tubes,"  Scien¬ 
tific  Report  No.  POLY- EE-80-007,  October,  1980). 

".V 

2)  Harmonic  Generation:  Project  is  completed.  A  technical 
report  is  issued: 

3)  Wave- Particle  Interaction  at  Cyclotron  Resonance 

The  motion  of  charged  particles  in  a  d.c.  magnetic  field  and  an 
rf  wave  is  governed  by  a  set  of  coupled  nonlinear  differential  equa¬ 
tions 

^  r(t)  =  v(t)  (3.1) 

3t  *(t)  *  §<£[?<t).t)]  +  x  tf[r(t).t]  +  £0[r(t)]})(3.2) 

which  are  not  subject  to  easy  solutions  in  general,  where  2(r,t)  and 
§(r,t)  are  the  oscillating  fields  of  the  EM  wave  and  S0(r)  is  the 
background  dc  magnetic  field. 

However,  the  usual  simplifications,  such  as  ignoring  the  per¬ 
turbation  in  the  trajectories  used  in  the  integration  of  the  equations 
of  motion  and  neglecting  the  oscillating  magnetic  field  of  the  rf  wave 
can  lead  to  significant  errors  at  cyclotron  resonances.  This  conclu¬ 
sion  is  manifested  in  the  previous  report  [Kuo  and  Cheo,  to  be  pub¬ 
lished]  which  used  a  self-consistent  resonant  trajectory  to  integrate 
the  equations  of  motion,  (3.1)  and  (3.2),  and  found  that  there  is 
significant  perpendicular  heating  and  no  parallel  heating  by  using 
ordinarily  polarized  wave  fields.  The  temporal  evolutions  of  the 


distribution  functions  and  perpendicular  temperature  are  also  obtained 
in  the  previous  works  for  a  simple  plasma  model,  i.e.,  uniform  back¬ 
ground  magnetic  field  and  normally  incident  EM  wave  (Kuo  and  Cheo, 

1980]  . 

In  the  presently  proposed  program,  electron  cyclotron  resonance 
heating  in  a  mirror  magnetic  field  will  be  investigated.  It  will  enable 
us  to  have  direct  comparison  between  our  theoretical  results  and  the 
experimental  observations  of  others.  The  bounding  effect  of  the 
trapped  electrons  in  the  mirror  field  on  the  resonance  heating  will  be 
investigated.  The  preliminary  results  [Kuo  and  Cheo,  1981]  have 
shown  that  the  trapping  effect  of  electrons  by  the  mirror  field  is 
indeed  important  to  the  resonance  heating.  It  is  because  the  phase 
variation  of  the  wave  seen  by  the  trapped  electrons  becomes  bounded 
due  to  the  finite  excursion,  all  the  trapped  electrons  can  receive 
resonant  interaction  with  the  heating  wave  at  a  frequency  equal  to 
either  the  fundamental  or  harmonics  of  the  mid-plane  cyclotron  fre¬ 
quency.  The  arising  of  the  additional  resonant  interaction  without 
experiencing  mis-match  between  wave  and  electrons  in  such  inhomo¬ 
geneous  magnetic  field  and  arbitrary  incident  angles  can  be  under¬ 
stood  from  the  fact  that  electrons  are  bouncing  back  and  forth  in  the 
mirror  field.  Wave-electron  interactions  are  thus  Doppler  shifted  by 
the  integer  multiples  of  the  bouncing  frequency  and  one  component 
always  contributes  to  the  resonance  interaction  without  mismatch. 
The  perpendicular  heating  rate  and  the  effect  heating  wave  upon  the 
confinement  of  the  electrons  will  be  determined.  It  is  also  our  plan 
to  derive  the  adiabatic  invariant  relations  of  motion  under  the  reso¬ 
nance  conditions  for  various  mode  type  of  heating  waves  [Kuo  et  al., 

1981] .  The  results  of  the  study  may  provide  information  about  the 
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saturation  level  of  the  heating.  Other  possible  mechanisms  such  as 
relativistic  effect  and  synchrotron  radiation  will  also  be  investigated. 

4)  Turbulent  Interaction  Between  Waves  and  Charged  Particles 

The  problem  of  turbulent  heating  of  plasma  by  laser  fields 
through  parametric  excitation  is  investigated.  The  main  objective  of 
this  investigation  is  to  understand  the  heating  profile  in  the  velocity 
space  of  electrons  in  order  to  predict  and  also  to  verify  the  experi¬ 
mental  observations ,  i .  e. ,  when  and  how  the  suprathermal  electrons 
will  be  produced  and  how  much  heating  energy  will  be  deposited  into 
the  main  body  of  the  plasma  electrons.  In  previous  analyses  [Kuo 
and  Cheo,  to  be  published],  the  ratio  Yb(0)/yt(0)  of  the  initial 
energy  deposition  rates  to  the  bulk  electrons  and  to  the  suprathermal 
tail  electrons  was  determined.  It  was  found  that  bulk  heating  can 
dominate  over  tail  heating  when  the  pump  field  is  sufficiently  strong 
E0  >  (4rai0Te)1/2. 

It  is  known  that  the  wave  particle  interactions  become  more 
effective  than  the  result  of  conventional  quasi-linear  diffusion  when 
the  turbulence  level  increases.  This  is  because  the  orbit  of  the 
charged  particle  is  constantly  perturbed  by  the  turbulent  wave  fields 
and  hence  the  resonant  interaction  between  waves  and  particles  are 
broadened.  If  the  turbulence  is  parametrically  excited  by  a  coherent 
pump,  as  considered  in  the  present  analysis,  the  particle  orbit  is  also 
perturbed  by  the  pump  field  in  addition  to  the  turbulent  fields.  The 
oscillation  of  the  particle  due  to  the  pump  field  will  also  broaden  the 
interaction  with  the  turbulent  fields.  This  is  due  to  the  fact  that  the 
phase  velocity  of  the  turbulent  waves  as  viewed  by  the  oscillating 
particles  is  Doppler  shifted  by  integral  multiples  of  ut0/k,  where  u»o  is 
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the  pump  frequency  and  k  is  the  wavenumber  of  the  turbulent  field. 
Therefore,  the  parametrically  excited  electron  plasma  wave  does 
interact  with  those  bulk  electrons  effectively.  Moreover,  the  produc¬ 
tion  of  electrons  at  several  phase  speeds  of  the  electron  plasma  waves 
is  also  predicted. 

The  results  of  the  previous  analysis  is  preliminary  in  the  sense 

\ 

that  only  the  ratio  of  initial  energy  deposition  rates  is  determined, 
and  also  the  case  of  E§/47tn0Te  »  1  so  that  the  possibility  of  the 
production  of  ultra-suprathermal  electron  to  the  anomalous  level  is  not 
considered. 

In  the  present  investigation,  temporal  evolution  of  the  ratio  of 
energy  deposition  rates  to  the  bulk  electrons  and  to  the  suprathermal 
and  ultra-suprathermal  tail  electrons  will  be  derived.  The  depen¬ 
dence  of  the  level  of  the  produced  ultra-suprathermal  electrons  on  the 
pump  power  will  be  determined  and  the  effects  of  the  inhomogeneity 
and  the  self-generated  dc  magnetic  field  in  the  background  plasma  on 
the  present  energy  coupling  mechanism  will  also  be  studied. 

5)  RF  Driven  DC  Current  Experiment 

Two  models  of  rf  driven  dc  currents  have  been  investigated. 
The  first  includes  collisional  effects  based  on  momentum  conservation. 
This*  technique  yields  the  following  dc  electron  and  ion  velocities 


,  k  e2 

ve  =  2  mfef  £2  xe 
e  e 

1  ku  e2 

vi  =  '  2  Sp"  E  Ti 


(5.1) 

(5.2) 
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The  resultant  dc  current  density  is  given  by 

-  (— Y 

c  \B„/ 


J  -  e_  nQ_ 


4V5t  k„ 


ut 


V5T  u  Jn 


3*  '  “ 

:  u>  u.J 


(5.3) 


The  second  model  is  based  on  retaining  the  second  order  vx  x  bx 
term  in  the  equations  of  motion. 

The  dc  current  density  in  this  case  is  given  by 

J  =  ejioju.  /B_  V  r  +  Si  (5 .4) 

4 -Jn  \B0  f  [  u)  ] 

Equations  (5.3)  and  (5.4)  are  valid  in  the  frequency  range, 

a  «  u»  <<  Qe  (5.5) 


The  equations  are  valid  for  electromagnetic  waves  propagating  in  the 
z-direction  where 


u»  =  u)(k)  (5.6) 

is  determined  by  the  dispersion  relation  of  the  specific  mode.  The 
details  of  the  above  results  were  presented  in  detail  in  the  previous 
proposal. 

Since  the  previous  proposal,  emphasis  has  been  placed  on  the 
completion  and  testing  of  experimental  facilities  and  the  design  of  the 
wave  launching  structure.  A  summary  of  the  experimental  program  is 
given  below. 

a)  Wave  Launching  Structure  for  Alf^en  Wave. 

The  problem  is  to  design  a  wave  launching  structure  that  will 
excite  an  Alftfen  wave  propagating  along  the  background  magnetic 
field  B0.  The  wave  will  be  coupled  magnetically  to  the  plasma  by 
means  of  an  artificial  transmission  line.  Magnetic  coupling  is  used 
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since  the  Alfven  wave  possesses  a  low  wave  impedance;  i.e.,  the 
wave  phase  speed  is  much  less  than  c. 

Basically  the  structure  will  be  constructed  as  shown  in  Figures 
(5.1)  and  (5.2). 

Figure  5.3  shows  an  equivalent  network'  for  the  artificial  trans¬ 
mission  line. 


Source 

Figure  5.3  -  Network  Representation  of  Artificial  Transmission  Line 

A  summary  of  pertinent  design  considerations  will  not  be  dis¬ 
cussed  . 

The  dispersion  relation  for  an  Alfven  wave  is  given  by 


(5.7) 


(5.8) 


where  B0  is  the  dc  background  magnetic  field  and  n0  is  the  ion 
density . 

We  see  from  (5.7)  that  for  a  given  background  magnetic  field 
and  gas,  that  the  only  experimental  variable  is  n0  and  u>.  These 
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parameters  must  be  chosen  with  some  care  so  as  to  make  the  experi¬ 
ment  feasible.  Considering  the  second  theoretical  model,  equation 
(5.4)  was  the  resultant  dc  current  density.  Inserting  Eq.  (5.7)  into 
Eq.  (5.4)  yields  the  following  expression  for  dc  current  density. 


T  -  e2  nft  1/2 /B_\2 
J  16tt  |j0  m.  \B0 ) 


(5.9) 


The  dc  current  density  can  be  maximized  by  making  n0  as  large  as 
experimentally  possible.  However,  for  a  fixed  u»,  the  wavelength  of 
the  Alfven  wave  becomes  smaller. 

It  is  not  desirable  to  have  the  wavelength  become  so  small  that 
the  construction  of  the  artificial  transmission  line  will  become  unrea¬ 
sonable. 

The  first  step  in  the  design  is  to  determine  the  values  of  the 

network  parameters  L,  C,  and  RT  .  R_  is  the  internal  impedance  of 

l,  g 

the  rf  source.  We  desire  to  simulate  a  slow  traveling  wave  on  the 
artificial  transmission  line  so  as  to  couple  quasistatically  to  the  Alfven 
wave.  The  coupling  is  obtained  at  the  locations  of  the  rf  coils.  The 
capacitors  serve  as  delay  elements  so  as  to  appropriately  adjust  the 
phase  of  the  signal  at  various  points  in  the  network.  If  the  phasing 
is  adjusted  so  that  the  time  phase  delay  At  between  adjacent  coils  is 
equal  to  the  time  it  takes  the  Alfven  wave  to  propagate  a  distance  Az 
we  will  have  a  maximum  coupling.  The  impedance  of  the  transmission 
line  is 


Z0  = 


(5.10) 
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To  eliminate  reflections  and  to  have  maximum  power  transfer  from  the 
high  power  rf  source  we  set 

zo  =  Rl  =  Rg  (5.1D 

This  prevents  a  standing  wave  from  being  set  up  on  the  transmission 

line.  Since  R  and  R,  are  fixed  by  the  available  equipment  Z0  is 
g  u 

known  and  an  equation  involving  L  and  C  is  given  by  (5.10).  The 
phase  velocity  on  a  transmission  line  is  given  by 


v 


P 


1 

/  L_  C_ 
,/  Az  AZ 


(5.12) 


v_  = 


=  f\=-^- 


p  ~  k 


(5.13) 


Az  is  relatively  arbitrary  except  that  it  must  be  chosen  as  a  realistic 
number  in  terms  of  the  experimental  setup.  Once  the  network  para¬ 
meters  L  and  C  are  determined  one  can  determine  the  currents 
through  the  industors  L  and  hence  the  generated  rf  magnetic  field  for 
a  given  input  power. 

The  next  step  in  the  design  makes  use  of  dispersion  relation 
(5.7)  given  a  relationship  between  tu,  \,  and  n0.  For  a  given  density 
the  value  of  the  rf  magnetic  field,  B  will  determine  the  dc  current 
density.  B  will  be  determined  by  the  rf  current  flowing  in  the  rf 
coupling  coils  of  a  given  geometry. 

b)  RF  System 

Presently,  two  relatively  high  power  rf  generators  will  be  avail¬ 
able  for  experiments.  The  first  was  described  in  a  previous  report. 
This  unit  is  capable  of  delivering  cw  rf  power  of  lkW  from  15MHz- 
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50MHz.  It  consists  of  basically  three  amplifiers.  The  first  stage  is  a 
low  power  class  A  solid  state  driving  amplifier  of  approximately  3W 
output.  The  second  stage  amplifier  is  a  class  B  solid  state  amplifier 
capable  of  delivering  140W  output  at  full  drive  power.  The  last  stage 
is  a  class  C  grounded  grid  amplifier,  tunable  from  10MHz-50MHz 
capable  of  delivering  1KW  CW  with  full  drive  power  of  approximately 
75W. 

The  second  high  power  rf  source  is  a  tunable  5KW  CW  oscillator 
operating  in  a  frequency  range  from  2MHz-18MHz.  It  contains  the 
appropriate  audio  drivers  and  modulators  to  operate  the  unit  as 
modulated  continuous  wave  (MCW),  if  desired.  It  can  also  be  oper¬ 
ated  CW. 

c)  Pulse  Modulator 

To  increase  the  ion  density  range  of  the  experiment,  a  high 
power  pulse  modulator  has  been  developed.  This  unit  was  described 
in  the  last  proposal.  The  unit  is  capable  of  delivering  peak  powers 
of  32  MW;  16KV  @  2KA  pulses.  The  average  power  capability  has 
been  upgraded  from  approximately  8KW  to  16KW. 

The  new  unit  can  operate  up  to  a  repetition  rate  of  500  Hz  with 
a  pulse  width  of  approximately  1  ps,  using  a  line  type  modulator  and 
a  hydrogen  thyratron  as  the  switch  tube. 

d)  Experiments 

After  the  artificial  transmission  line  is  installed,  experiments  will 
proceed  in  the  following  manner. 

The  plasma  will  be  generated  by  ECRH  using  1KW  CW  microwaves 
at  2.45GHz.  The  rf  units  and  plasma  parameters  will  be  adjusted  to 
produce  a  dc  current  flow.  The  steady  state  background  density  n0 
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is  relatively  low  and  hence  limits  the  maximum  allowable  dc  current 
flow. 

The  S  band  microwaves  cannot  produce  a  high  density  plasma 
since  the  plasma  tends  to  screen  out  additional  microwave  power  when 
the  plasma  density  increases  to  a  value  where  the  corresponding 
plasma  frequency  is  equal  to  the  source  frequency  of  2.45GHz.  This 

_3 

critical  density  corresponds  to  7.4  x  10 10  cm  .  This  density  appears 
in  the  vicinity  of  the  microwave  coupler,  and  in  fact  the  density  is 

-3 

significantly  lower  in  the  test  region;  109  cm  .  To  investigate 
higher  dc  currents  a  method  to  enhance  the  density  must  be  devel¬ 
oped.  To  do  this  the  32MW  pulse  modulator  will  be  used  to  enhance 
the  density.  The  method  used  to  carry  out  this  experiment  is  as 
follows.  A  steady  state  plasma  is  produced  using  ECRH  with  the 
S-band  microwave  system.  A  pulsed  longitudinal  discharge  is  initiated 
with  the  pulse  modulator.  This  ionizing  high  current  pulse  produces 
a  high  density  plasma  afterglow  which  decays  to  its  steady  state 
density  in  a  few  collision  times. 

If  the  rf  drive  is  on  during  the  entire  experiment,  a  dc  pulse 
current  due  to  the  rf  drive  will  be  produced  at  a  point  where  the 
density  n0  is  such  that  the  Alfven  wave  dispersion  relation  is  satis¬ 
fied.  See  Figures  (5.4)-(5.7). 
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